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Unlike bones, behavior does not fossilize, so it is hard to infer the evolutionary history of social traits. However, we have 
shown elsewhere that Bayesian phylogenetic methods allow the investigation of ancestral states and models of evolution 
of social grouping behavior in primates. Here, we extend this analysis to another significant aspect of primate social life, 
which may be subject to different evolutionary pressures - mating systems. We show that mating systems evolved from a 
polygynandrous state at the root of the phylogeny to the two derived states of harem-polygyny and monogamy. Unlike 
social organization, where there were no transitions from uni-male groups to pairs, here we found positive transition rates 
from both polygynous mating states into monogamy. There were no transitions out of monogamy to another mating 
state. Both derived mating systems evolved late in primate evolution. Nocturnal primates remained solitary foragers 
while their mating systems evolved from polygynandry to harem-polygyny and monogamy. However, among diurnal 
primates the derived mating states evolved at the same time as the derived states of social organization. 



Although rare among mammals, sociality is common among pri- 
mates. 1 In a recent study we showed how group living evolved 
in primates following a switch to diurnal activity patterns. 2 
However, primate societies are shaped not only by the groups 
that animals live in, but also by the mating strategies that they 
employ. Therefore, for a full understanding of the evolutionary 
history of primate social systems both need to be investigated. 

We have shown previously that for social traits with a strong 
phylogenetic signal, such as social organization or dispersal pat- 
terns, it is possible to infer models of evolution, transition rates 
between states and ancestral states at nodes on the primate phy- 
logeny. 2 Here we use Bayesian phylogenetic methods 3 to infer 
ancestral states and evolutionary pathways for mating systems 
in primates, complementing our previous investigation of social 
organization and allowing the comparison of these two features 
of primate life. 

An analysis using Ape 4 in R gave a Lambda value of 0.996, 
not significantly different from 1 (p = 0.679), suggesting a strong 
phylogenetic signal in the mating system data among primates. 

We ran RJ MCMC analyses in BayesTraits 5 with mating sys- 
tem classified as both a binary trait (polygyny/monogamy) and 
a three state trait (harem-polygyny/polygynandry/monogamy). 
The analysis for the binary classification supports polygyny as 
the ancestral primate mating system at the root of the phylogeny 
(polygyny mean probability = 0.974 +/- 0.001; monogamy mean 
= 0.026 +/- 0.001). For the three state trait the root was polygy- 
nandry (polygynandry mean = 0.853 +/- 0.003; harem-polygyny 
mean = 0.109 +/- 0.002; monogamy mean = 0.038 +/- 0.001). 
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The best fitting model from the RJ procedure (with 71% of 
the posterior probability distribution) revealed that from polygy- 
nous mating at the root of the phylogeny there was a strong tran- 
sition rate to monogamy, but a zero reverse rate (Fig. 1A). We 
also split the polygynous mating state into harem-polygyny and 
polygynandry to understand the dynamics of mating change bet- 
ter. 6 From polygynandry at the root of the phylogeny there were 
strong transition rates into both harem-polygyny and monog- 
amy, and a non-zero but weaker rate from harem-polygyny to 
monogamy (best fitting model from the RJ procedure, with 33% 
of the posterior probability distribution). All other rates were zero 
(Fig. IB). 

Across the whole posterior probability distribution, for the 
three state analysis, the mode transition rate out of monogamy 
was zero, while the transitions from harem-polygyny to monog- 
amy and back to polygynandrous mating were zero for 18% and 
41% of the time respectively (Fig. IB). Transitions from polygy- 
nandrous mating to the other two states were never zero. 

Harem-polygyny evolved earliest among strepsirrhines, first 
at the root of the Loris sub-family Perodicticinae (~42mya) and 
later at the root of the Lepilemurs (-36mya) (Fig. 2). Monogamy 
evolved in a number of Lemur families from -28mya. Among 
haplorrhines monogamy emerged first at the root of Callicebus 
(-26mya), followed by Aotus and the Callitrichids (~22mya) 
and then Gibbons (~19mya). Harem-polygyny evolved later 
among Anthropoids, at the root of the Trachypithecus genus 
(-I6mya), and later still at the root of the Cercopithecus genus 
(-llmya). 
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Figure 1. Model of evolution of primate mating systems showing posterior distribution of transition rates between states. A) Monogamy and po- 
lygyny. B) Monogamy, harem-polygyny and polygynandry. Thickness of arrows reflects proportion of time the transition rate is not zero. Z denotes a 
zero transition rate as a proportion of posterior probability distribution. A dashed line denotes a zero transition rate in the RJ derived model. Graphs 
show posterior probability distribution of each transition rate. 



Like social organization 2 mating system data in primates show 
a strong phylogenetic signal, indicating that it is possible to make 
inferences about the evolution of this trait across the primate tree. 
This also suggests that history plays a significant role in the cur- 
rent distribution of both these social traits across species. The 
model of evolution for primate mating systems is similar to that 
for primate social organization, 2 once multi-male/multi-female 
groups were established, with transitions to both monogamy and 



harem-polygyny. The pathways differ in that monogamy was 
found to evolve from both harem-polygyny and polygynandry. 

As with social organization, once monogamy was established 
(or pair living, in the case of social organization) there were no 
further transitions. This may be because the cognitive changes to 
enable the behavioral co-ordination required for stable monogamy 
are hard to reverse, 7 and also because the factors leading to the 
evolution of monogamy persisted over time. 8 ''' Harem-polygyny 
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Figure 2. Primate phylogeny with ancestral states for mating systems derived from RJ MCMC model of evolution. The tree topology is the maximum 
clade credibility tree from the lOkTrees Project 13 posterior distribution with branch length drawn proportional to time. Branches and tips are colored 
for polygynandry (red), harem-polygyny (orange) and monogamy (pink) where the combined probability of the state and the branch is greater than or 
equal to 0.7. Where the combined probability is less than 0.7 the branch is gray. 
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also evolved directly from polygynandry, but there were back 
transitions. In contrast, there were no direct transitions into 
harem-polygyny from monogamy suggesting that rather than a 
"combination of pairs" 10 harem-polygyny can be better described 
as "degraded" polygynandry. An interesting question for future 
work is the extent to which this is due to ecological and physi- 
ological conditions, such as small female groups and asynchro- 
nous estrus, enabling a single male to monopolise mating within 
a group of females. 1112 Co-evolutionary models could be used to 
test the factors leading to the emergence of harem-polygyny. 

The derived states of mating system and social organiza- 
tion evolved at a similar time across diurnal primate clades such 
that they matched through most of primate evolutionary his- 
tory. However, among those primates that remained nocturnal, 
solitary social organization persisted even when mating systems 
changed from the ancestral state of polygynandry to the derived 
states of harem-polygyny or monogamy. 

Methods 

Model testing was performed across a Bayesian posterior distribu- 
tion of 10,000 ultrametric primate trees derived from genetic data 
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